On the formulation of functional theory for pairing with particle number restoration 
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The restoration of particle number within Energy Density Functional theory is analyzed. It is 
shown that the standard method based on configuration mixing leads to a functional of both the 
projected and non-projected densities. As an alternative that might be advantageous for mass 
models, nuclear dynamics and thermodynamics, we propose to formulate the functional in terms 
directly of the one-body and two-body density matrices of the state with good particle number. Our 
approach does not contain the pathologies recently observed when restoring the particle number in 
an Energy Density Functional framework based on transition density matrices and can eventually 
be applied with functionals having arbitrary density dependencies. 
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I. INTRODUCTION 

Energy Density Functional (EDF) methods provide a 
universal framework to describe nuclear structure, nu- 
clear dynamics or thermodynamics. Tremendous ad- 
vances have been made in the last few decades on the 
practical application of EDF methods to nuclei [1] . Still, 
despite their long success, some of the fundamental as- 
sumptions made to justify the usual strategies how the 
EDF techniques arc constructed and used for nuclear sys- 
tems have not yet been satisfactorily clarified. Most, 
if not all, EDF approaches break as many symmetries 
of the nuclear Hamiltonian as possible: translational, 
rotational and U(l) symmetry in gauge space, among 
the most important ones. In fact, the exploitation of 
symmetry breaking in nuclei is strongly motivated by 
experimental observations. For instance, the appear- 
ance of highly collective rotational bands in spectroscopic 
data clearly points to the existence of deformed intrin- 
sic states in many nuclei [2]. Similarly, there is evidence 
that pairing can be often treated by explicitly breaking 
the U(l) gauge symmetry of eigenstates of the particle- 
number operator, like for instance in a Bardeen-Cooper- 
Schrieffcr (BCS) or Hartrec-Fock-Bogoliubov (HFB) ap- 
proach [2, 3]. Nuclei are, however, finite systems and 
methods like BCS or HFB do not properly treat quan- 
tum fluctuations of the order parameter associated with 
the broken symmetry [3]. These fluctuations can be in- 
corporated either by a statistical treatment of the order 
parameter, or by the restoration of the relevant symmetry 
[3]. The concept of symmetry breaking and restoration 
stands out as the tool of choice within the EDF frame- 
work. 



It has, however, been recently shown that restoration 
of symmetries has to be handled with great care in an 
EDF framework [4-7]. In particular, the configuration 
mixing within a Multi-Reference (MR) EDF approach 
might lead to serious practical difficulties that can, how- 
ever, eventually be cured [6, 8]. Besides compromis- 
ing applications, these difficulties have clearly pointed 
out the necessity to clarify the theoretical framework on 
which the theory can be build. 
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The discussion in the present paper is restricted to 
ground-state properties and to particle-number projec- 
tion, for which detailed analyses have been recently 
made. This case is the simplest situation where patholo- 
gies of the MR-EDF approach have been observed [4] , an- 
alyzed and regularized [6, 8]. The first goal of the present 
work is to provide an alternative analysis of the EDF 
theory using configuration mixing to restore symmetries 
without and with the regularization. It will be shown 
that neither the non-regularized nor the regularized func- 
tionals can straightforwardly be interpreted in terms of 
the densities of projected or non-projected states. Start- 
ing from this analysis, the second intent of this work is 
to propose an alternative way to introduce a functional 
theory that is U(l) symmetry conserving, and that with- 
out making use of the Multi-Reference technique. Our 
approach remains close to the Hohcnbcrg-Kohn [9] and 
Kohn-Sham [10] framework and uses a projected state as 
an intermediate trial state to construct the components 
of the functional. This approach avoids the difficulties 
recently encountered in MR-EDF approaches and can be 
applied also with functionals that cannot be safely em- 
ployed within the standard MR-EDF approach, as for 
example functionals with non-analytical density depen- 
dences. 
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II. PARTICLE NUMBER RESTORATION 
WITHIN EDF THEORY: STANDARD 
APPROACH 

The strategy to obtain a functional for pairing includ- 
ing particle-number restoration has been extensively an- 
alyzed recently [5, 6, 8, 11] and we only give here the 
main steps necessary for our discussion. Following these 
references, in this section we will consider a specific class 
of functional form that will be sufficient for the present 
discussion. 1 At the so-called Single- Reference (SR) level, 
a quasi-particle (QP) vacuum of Bogoliubov type |<£>o) 
is used to construct the normal and anomalous density 
matrices, denoted by p and k, that serve to construct the 
functional. The energy is then written as 



£sr[$o] = £sr[p,k,k*} 
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where v pp and v KK denote the effective vertices in the 
particle-hole and particle-particle channels. Standard 
SR-EDF can be schematically seen as the sequence 



(p,K,K*) 



£,SR ■ 



(2) 



The price to be paid for incorporating pairing with a 
rather simple functional is to use an intermediate state 
$o) that is not an eigenstate of particle number. In a 
second step, the symmetry can be restored projecting out 
the component with N particles 



Itfjv) = P N \$o), 



(3) 



where P N denotes the particle number projection oper- 
ator defined through [2, 3] 
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The expectation value of any operator O that conserves 
particle number can then be expressed as 



(* N \Q\V N ) _ 

(^n\^n) 
where the shorthand 
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(5) 



27T {^> N \^ N ) 



(G) 



1 Note that, none of the currently used SR-EDF functionals be- 
longs to this class as they have non-analytical density depen- 
dences. The form (1) is the only one (restricting ourselves here 
to bilinear functionals) for which the recently proposed regular- 
ization applies. 



has been introduced. Here ip denotes the gauge angle, 
whereas |<& v } = e lipN |<& ) refers to the state |$ ) rotated 
in gauge space. The kernel entering in the integral of 
Eq. (5) corresponds to the transition matrix element of an 
operator between two quasi-particle vacua. One can then 
take advantage of the Generalized Wick Theorem (GWT) 
[2, 12] to express the kernel in terms of the transition 
density matrices 
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For instance, when O is a two-body Hamiltonian, the 
two-body interaction v entering in (5) takes a form sim- 



ilar to Eq. (1) with v 
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and where the den- 



sities, p and k are replaced by the corresponding tran- 
sition densities, Eqns. (7-9). Guided by the Hamilto- 
nian case, the energy functional associated with particle- 
number restoration is usually defined through 



£ N [V N } = dtpS, 
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This energy functional is a special case of a so-called 
Multi- Reference EDF (MR-EDF). The present strategy 
to restore symmetries in an EDF framework deserves 
some further remarks: 

• The expression (10) is postulated having in mind 
the Hamiltonian case. However, the MR-EDF the- 
ory should not be confounded with the expecta- 
tion value of a Hamilton operator. In particular, 
an energy functional has much more flexibility re- 
garding the functional form of the energy kernels in 
Eq. (10), which can be used for the efficient mod- 
eling of in-medium correlations. 

• The construction of the MR-EDF, Eq. (10), from 
the SR-EDF by simply replacing the normal and 
anomalous density matrices in the SR EDF by the 
corresponding transition density matrices is postu- 
lated by analogy to the GWT. While it appears 
rather natural, it was shown recently that this 
strategy to construct the MR-EDF might lead to 
an ill-defined functional that exhibits divergencies 
and jumps in practical applications [5, 6, 8]. While 
a solution to this problem has been proposed and 
applied in Refs. [6, 8], a consistent framework for 
MR-EDF approaches is still missing. A clear illus- 
tration of this is the ongoing debate about which 
densities should enter in the functional [13], as well 
as the recently recognized impossibility to use non- 
integer powers of the transition density in MR en- 
ergy functionals [11]. 
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• The very notion of symmetry restoration within an 
EDF framework remains to be clarified. For in- 
stance, it has been shown recently [7] that also reg- 
ularized MR energy functionals may in general not 
transform as an irreducible representation of the re- 
stored symmetry, unless additional constraints are 
introduced. 

In the present section, we will further analyze the way 
the MR-EDF is constructed and the possible sources of 
difficulties. For simplicity, we restrict ourselves to a case 
where the two-body effective interaction kernels entering 
Eq. (1) are independent of the densities. 

A peculiarity of particle-number projection is that the 
canonical basis of the original state |<f>o) an d of the ro- 
tated states \<&u>) are the same when making a suit- 
able choice of the Bogoliubov transformation between 
quasi-particle operators of these states. Accordingly, the 
canonical base of the projected state \^>n) is also the 
same as the one of the original reference state | <&()}• hi 
the following, we will implicitly assume that densities are 
written in this canonical basis in which we have 



Pij = 5 H n i . 



whereas the energy £n takes the form 

,■•2 71 



V 0* 



(11) 



£n[*n] = V t u / dip n? v Af N (0, if) 
Jo 

X k£\°£jV n (0,ip).(12) 



After a lengthy, but straightforward calculation, the en- 
ergy functional can be expressed as 



£n[*n] = J2 U 
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where n f are the occupation numbers: 
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and Rfj k[ corresponds to the two-body density matrix 



(15) 



of the projected state. They can be expressed in terms 
of the gauge angle integrals as 



,N 
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and 
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(17) 



Equation (13) is rather enlightening with respect to the 
physical content of present MR-EDF calculations. In- 
deed, if one neglects the last two terms in Eq. (13), one 
sees that the functional associated with the projected 
state can be regarded as a functional of the one- and 
two-body components of this very state. Similarly, if one 
uses the same effective interaction v pp = v KK , then the 
last two terms of Eq. (13) recombine and the two-body 
component can be recognized, thanks to the relation 



dip (n. 



<°^-f)A6v(0^). (18) 



However, when using different effective vertices v pp ^ v KK 
in the particle-hole and particle-particle channels, or 
when using vertices v pp or v KK that cannot be written 
as an antisymmetrized matrix elements of the two-body 
force, then the identification of the energy as a functional 
of one- and two-body density matrices of the projected 
state cannot be made anymore. Instead, it can only be 
written as a functional of the transition density matri- 
ces. 2 This subtlety is intimately connected to the pres- 
ence of pathologies encountered in MR-EDF calculations. 
Indeed, the last two terms in Eq. (13) are nothing but the 
ones at the heart of the difficulties to construct a well- 
defined MR-EDF theory. As discussed in Refs. [5, 8], 
for near-orthogonal states ($o|^V) — there is at least 
one n® v and the corresponding nf® and K^f that all go 
to infinity. As a consequence, the two terms can sepa- 
rately become larger than any physical scale in the nu- 
cleus. They do, however, recombine to a well-behaved ex- 
pression when a Hamiltonian is used, i.e. when v pp = v KK . 
Without taking specific care of these terms in the restora- 
tion of symmetry within the functional framework, there 
is a spurious contribution that leads to discontinuities 
and divergences when plotting the particle-number pro- 
jected energy as a function of a collective coordinate. 



2 We recall that the expectation value of the two-body operator 
in a projected state can be written as a functional of the two- 
body density of this state, or, fully equivalently, as a functional 
of the one-body density matrices. This property does not hold 
for general functionals that are constructed without reference to 
an underlying Hamiltonian. 
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A. MR-EDF theory with regularization 

A strategy to construct a well-behaved MR-EDF the- 
ory proposed in Rcfs. [6, 8] is to remove terms that might 
not properly rccombinc in the MR-EDF approach in such 
a way that the spurious contamination is removed with- 
out touching the physical content of the functional. The 
resulting functional then takes the form (technical details 
are given in appendix A) 



£n[^n] 



Z 

- V v pp - R r 

4 2-^< "33 3 

i 

lJ2 v ^K( l - n ?) + Sn ^} > ( 19 ) 
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where Sri; 



n® is the difference between the oc- 



cupation number of the level i in the projected and the 
non-projected state. 

Expression (f 9) is of particular interest for the follow- 
ing discussion regarding the construction of energy func- 
tional theory. First, let us remark that, compared to the 
previous form (13), the gauge space integrals are now 
hidden in the components of the one- and two-body den- 
sity matrices of the projected state. In addition, the last 
two lines of Eq. (19) are also functionals of the occu- 
pation numbers in the original non-projected state. 
The analysis of the regularization procedure to remove 
spurious contribution to the MR-EDF method [6, 8, 11] 
suggests that these terms will always be well-behaved. 

An example for a deformation energy curve obtained 
from a particle-number projected MR-EDF calculation 
with the Skyrme interaction SIII and a pairing functional 
of volume type is shown in Fig. 1 (dashed line). The MR- 
EDF is numerically calculated using expression (10) and 
the Fomenko discretization procedure of the gauge-space 
integrals described, for instance, in Ref. [8]. Here, 199 
discretization points have been used. This large number 
is necessary to resolve the discontinuities that stem from 
the spurious contribution to the non-regularized MR- 
EDF [8]. As in Ref. [8], the Lipkin-Nogami procedure 
is used in the minimization of the energy of the state 
|$o)- The solid line corresponds to the MR-EDF method 
with the regularization proposed in [6]. In this Figure, 
we also show the results (filled circles) obtained using 
directly the expression (19) that has been proven above 
to be analytically equivalent to the regularized MR-EDF 
functional. Note that, in the latter case, we have used 
a method called hereafter "recurrence method" to com- 
pute the projected occupation numbers and components 
of the projected two-body densities. This method is de- 




FIG. 1. (Color online) Particle-number restored deformation 
energy curve of ls O as a function of quadrupole deformation 
/?2 and calculated within standard MR-EDF technique using 
Projection After Variation (PAV) with SIII and a delta pair- 
ing interaction before (blue dashed curve) and after regular- 
ization (black solid curve). The red filled circles correspond 
to the result obtained using directly Eq. (19) (see text). To 
compare with previous work [8], the Coulomb exchange con- 
tribution has been subtracted from the energy. 



scribed in detail in appendix B. Although the use of gauge 
angle integration would have given exactly the same re- 
sults, this method has the advantage to be very simple, 
numerically efficient and to not make use of transition 
density matrices. As expected, the energy obtained with 
expression (19) exactly matches the one obtained using 
the regularized MR-EDF functional. This formulation 
provides a new and alternative insight into the content 
of particle-number restored energy functionals. 



B. Critical analyses of standard method 

As discussed above, specific regularizations in MR- 
EDF functionals are needed to avoid discontinuities such 
as the jumps appearing in Fig. 1. At this point, even 
with the regularization, two important problems remain: 

(i) Terms that have a non-analytical density depen- 
dence, for example a non-integer power such as of- 
ten used in parameterizations of the Skyrme and 
Gogny interactions, cannot be regularized with the 
procedure proposed in Ref. [6]. Indeed, the func- 
tional itself becomes in that case multivalued in 
the complex plane and cannot be properly de- 
fined [5, 11]. 

(ii) A second issue illustrated in Eq. (19) is that the 
last two terms are not only a functional of the oc- 
cupation numbers of the projected state, but also 
of the occupation probabilities of the original refer- 
ence state | } - Accordingly, the energy remains a 



5 



functional of the density of the quasi-particle vac- 
uum that is not an eigenstate of particle number. 
This raises the question which density, i.e. pro- 
jected, transition, or non-projected can be used to 
construct a functional for MR calculations. 

In the following, we show that both (i) and (ii) can 
eventually be avoided by changing the strategy to con- 
struct the functional for pairing that accounts for particle 
number restoration. 



III. DISCUSSION ON EDF THEORY FOR 
PAIRING WITH PARTICLE NUMBER 
RESTORATION 

Let us now discuss the critique (ii) made above con- 
cerning the components of the projected energy func- 
tional. In most functional approaches, an intermediate 
state is introduced to construct densities that are used 
to minimize the energy. This is the case in usual DFT 
or at the SR-EDF level where the trial state is a Slater 
determinant or a quasi-particle vacuum. When restoring 
the symmetry in a MR-EDF framework, then, according 
to Eq. (19), the projected state can be almost regarded 
as an intermediate many-body state from which the one- 
and two-body density matrices used to define the func- 
tional are obtained. 

However, due to the presence of in the energy, this 
functional happens to depend on components not only 
of the projected state, but also of the original reference 
state. A slight modification, however, can easily restore 
the unique dependence of the functional on the projected 
state. If, for instance, the following replacements 

{n^nf — SriiSrii) => n^nf 

(20) 

[nf(l-nf) + =► nf(l-nf), 

are made in Eq. (19), then the strategy of standard DFT 
to construct the EDF as a functional of a density of an 
auxiliary state, the projected state here, is recovered. 3 

The use of a projected product state a auxiliary state 
has the advantage that it allows to treat pairing in a 
particle-number conserving framework. An illustration 
of a result obtained taking into account this modification 
in Eq. (19) is shown in Fig. 2 and compared to the orig- 
inal curve. This figure illustrates that the small change 
in the functional does not significantly modify the en- 
ergy landscape. This is indeed not unexpected since the 
difference drii (resp. SriiSrii) is likely to be much smaller 
than (resp. nf nf). 



3 This does not mean, however, that we recover a theory that is 
equivalent to DFT. Indeed, at this stage, the functional (19) is 
still a functional of the two-body density matrix. However, as 
will be discussed below, for the specific case of particle-number 
projection, the two-body density matrix is itself a functional of 
the one-body density matrix. 
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FIG. 2. (Color online) Particle-number restored deformation 
energy surface of 18 calculated using Eq. (19). The dotted 
curve is obtained by making the replacement (20) in Eq. (19). 

By making the simple modification (20), the EDF 
framework can now be interpreted as a functional of the 
projected-state densities. Indeed, the state with good 
particle number can now be regarded as the auxiliary 
many-body state that provides the quantities used to 
construct the EDF. Similarly to Eq. (2), the correspond- 
ing theory can be regarded as a new sequence 

* N => (p n ,R n ) => £n- (21) 

We note in passing that the slight modification (20) does 
not break the shift invariance of the energy functional 
discussed in Refs. [5, 8]. At this point, let us make a 
few further important remarks: 

• The functional form (19) is not completely surpris- 
ing. Indeed, our starting point, Eq. (1), is very 
close to a form one would have obtained by taking 
the expectation value of a two-body Hamiltonian. 
In the case of an energy functional calculated as 
the expectation value of a genuine Hamiltonian op- 
erator, the energy can be written as a functional of 
one-body transition densities, but also as a func- 
tional of projected one- and two-body (and eventu- 
ally higher-order) densities. This freedom is lost in 
the functional framework, where a choice has to be 
made for cither one or the other. MR-EDF follows 
the former strategy, whereas the latter has not been 
explored yet. For a regularized bilinear functional, 
the differences between both formulations remain 
very small, see Fig. 2. 

• Expression (19) contains not only one-body but 
also two-body matrix elements and might appear 
out of the scope of a density functional approach 
aiming at replacing the original A-body problem 
by a functional of the one-body density matrix 
only. Indeed, in the Hamiltonian case, the expec- 
tation value of any two-body Hamiltonian for any 
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state can directly be regarded as a functional of the 
one- and two-body densities of this state. Density 
functional theories are introduced to avoid the ex- 
plicit use of two-body and higher density matrices. 
Therefore, by itself, the use of a functional of the 
two-body density might appear useless. The im- 
portant simplification here stems from the fact that 
these densities are constructed from a very specific 
class of states, namely projected product states. 
For instance, we have shown recently that the two- 
body density matrix elements can eventually be 
written as an explicit functional of the one-body 
density under some approximation [14]. Accord- 
ingly, while two-body density matrix elements are 
used to get a compact expression in Eq. (19), this 
functional can truly be considered as a functional 
of the projected state one-body density consistently 
with density matrix functional theory, such that the 
sequence become 

*w P N =>• £n- (22) 

• When making the replacement (20) in Eq. (19), 
then the functional directly incorporates symmetry 
breaking and its restoration in a single step, con- 
trary to standard approaches in EDF theory. From 
that point of view, it could be seen as a " Symmetry- 
Conserving" EDF theory. 4 

• It is quite interesting to note that the MR-EDF 
can already almost be regarded as a functional of 
the components of the projected state. While this 
was hidden in formula (10), it becomes evident in 
Eq. (19). In particular, as noted in Refs. [6, 8, 11], 
there exists some flexibility in the rcgularization of 
the pathologies of the MR-EDF. It is possible to 
slightly modify the original prescription (A4-A5), 
such that the rcgularization automatically leads 
to (20). In that case, the method based on the 
use of MR-EDF and "symmetry conserved" EDF 
framework are strictly equivalent. As an important 
consequence, while the use of techniques inspired 
from configuration mixing was unclear within a 
functional framework, we give here evidence that 
it can be formulated consistently in a functional 
framework. It is, however, worth mentioning that 
while this connection can be made only in the sim- 
ple functional form given in Eq. (1), most func- 
tionals currently used do not allow their controlled 
usage in an MR EDF framework. 

• Finally, it is important to mention that this equiv- 
alence holds only true for the schematic bilinear 



4 It should be, however, kept in mind that the present functional 
only takes care of the restoration of U(l) gauge symmetry while 
others still remain broken. 



functional given by Eq. (1) with two-body vertices 
independent on the density. If density dependent 
terms are present in the functional, like in all cur- 
rently used parameterizations of the EDF, such an 
equivalence does not exist anymore. Note, however, 
that, in this case, a safely usable MR-EDF cannot 
be constructed anymore due to the absence of a 
rcgularization scheme. In Eq. (19), one then ob- 
tains a functional that remains closer to the spirit of 
DFT based on the Hohenberg-Kohn theorem than 
the usual MR-EDF approach. Indeed, in the HK- 
theorem-based DFT, the functional is constructed 
from the density matrices of the correlated (i.e. in 
our case projected) state. As we will illustrate be- 
low, on the contrary, the alternative formulation 
proposed here that treats both symmetry breaking 
and restoration simultaneously can still be applied 
for functionals that cannot be regularized in a MR- 
EDF framework. 



A. Constraints on the symmetry-conserving 
functional 

If the standard projection method is used as guidance 
to construct the functional, then the form of the func- 
tional is almost entirely constrained. Indeed, this corre- 
sponds to use Eq. (13) or eventually Eq. (19) as a starting 
point. Eq. (20) corresponds to a specific choice. Here, 
we discuss whether alternative choices can be made for 
the last two lines of Eq. (19). At present, it is not clear 
if, within the functional framework, a unique prescrip- 
tion of the functional form exists. Nevertheless, one can 
propose a few rules to better constrain its form. Let us 
assume a more general prescription than Eq. (20) 

(nfnf-SmSn,) => Fg, 

(23) 

{nf{l-nf) + 5 ni 8m) =► G% , 

where F N and G N are the unknown quantities. 
Let us specify some rules to constrain them: 

o Sum-rule: When v pp = v KK , then the last two 
terms in Eq. (13) should recombine to give R~ = 
nf . Accordingly, it seems reasonable to impose 

Fg + G»=n?. (24) 

o No-pairing limit: Slater determinants belong 
to the Hilbert space spanned by projected states. 
Consequently, one can interpret the functional for 
particle-number projected wave functions as a gen- 
eralization of the SR-EDF theory expressed for 
Slater determinant, i.e. 

SnI^n] =>■ £s*[*sd] . (25) 

as $at — > $sd- ^sd denotes any Slater determi- 
nant. As a consequence, in this limit, we should 
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have 



N 



JV 



0. 



(26) 



o Large TV limit: In the limit of infinite particle 
number, the projected state and the reference state 
should become identical (for instance Sn® ==> 0). 
Accordingly, we do expect 



lim 



N- 



F N - n N n N 



lim 



N- 



Cx, v 



<(1 



(27) 



These three constraints significantly reduce the freedom 
of choosing the form of the functional that can be used. 
The prescription (20) naturally fulfills all of them. 



B. Can we use terms with non integer power of the 
density? 

When the effective two-body vertex depends explic- 
itly on the density, then the energy cannot be directly 
mapped on Eq. (13). If the density dependence is in inte- 
ger powers of the density, then one could eventually gen- 
eralize the derivation of Eq. (13) to three-body or even 
higher-order effective interactions. For all other forms 
of the density dependence, such as the widely used non- 
integer powers of the density, there is no way to deduce an 
equivalent expression because the integration over gauge 
angles cannot be uniquely defined from a mathematical 
point of view [5, 11]. It is worth to mention that the same 
difficulty appears when the Coulomb exchange term is 
approximated using the Slater prescription. Above, we 
have shown that, with a slight change in the functional 
used in the standard MR-EDF method, one obtain an 
EDF that can can be interpreted consistently within the 
usual functional approach where the projected state be- 
comes a trial wave function to construct the ingredients 
of the functional. 

Guided by the setup of functional (13), the most nat- 
ural and simple way to extend the SR-EDF functional 
using density dependent two-body effective vertices with 
non-integer powers of the density is to directly replace 
the density entering in the effective vertex by the density 
of the projected state, i.e. 



v pp [p] => v pp [p N ] 



■'[f>] 



-kk r „n 



V 



(28) 



Again, by doing this, we ensure that the functional used 
for the projected state is consistent with the one used in 
the no-pairing case (Eq. (25)) and in the large- N limit. 

In Fig. 3, the deformation energy curve obtained by 
using Eq. (28) is compared to the result deduced from 
the standard non-regularized MR-EDF procedure using 
Eq. (10). The SLy4 effective interaction used here con- 
tains density dependent terms with non-integer powers 
i.e. p 2+1 / 6 . Note that in this case, the MR-EDF cannot 




FIG. 3. (Color online) Same as Figure 1 when the SLy4 ef- 
fective interaction is used in the particle-hole channel. The 
dashed line corresponds to the non-regularized MR-EDF re- 
sult directly obtained by gauge angle integration using 199 
points in the discretization. The filled circles correspond to 
the result obtained using the Symmetry-Conserved EDF us- 
ing Eq. (19) and the prescription (20). 



be regularized. The new alternative method we propose 
here, however, does lead to a perfectly well behaved en- 
ergy curve. 

In our approach, the main difference between effective 
interactions that depend on non-integer powers of the 
density and those depending only on integer powers of 
the density, is that while in the latter case one might 
eventually recognize terms coming from three-body or 
four-body forces and so on, this is impossible in the for- 
mer case. It should, however, be kept in mind that the 
use of effective interactions should be regarded more as 
a guidance for the set-up of the actual form of the func- 
tional, and not as a prerequisite for the functional ap- 
proaches as such. 

The example of non-integer powers of the density 
shows that functional theory including pairing and parti- 
cle number restoration and extending the usual SR-EDF 
approach, but without using the MR-EDF framework, 
can eventually be defined for rather general class of func- 
tionals if the strategy to construct the functional pro- 
posed here is followed. Let us add a few remarks: 

• It is important to realize that for particle-number 
projection the present strategy becomes equiva- 
lent to the MR-EDF one when the regularization 
is slightly modified compared to the one originally 
proposed in Ref. [6], i. e. the present strategy and 
the modified regularized MR EDF calculation will 
give the same energy for regularizablc functionals. 
For those, it should therefore be more regarded as 
an alternative way of implementing MR-EDF ap- 
proach to particle-number projection than as a new 
framework. 

• With the present strategy, one will never have prac- 
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tical difficulties in applying the method to rather 
general and complex forms of functionals. How- 
ever, some effort has been made recently to out- 
line the constraints that a functional should ful- 
fill to be truly regarded as a symmetry-conserving 
functional [7]. While these constraints are even 
partially unknown, one might anticipate that they 
will significantly restrict the functional form that 
might be used. We are therefore facing the fol- 
lowing dilemma: from condensed matter physics, 
we know that the powerfulness of DFT comes from 
the large flexibility in choosing the functional form. 
Putting too many formal constraints will make it 
increasingly difficult to model the relevant physics 
with a tractable functional. In particular, one can 
already see from Ref. [7] that a functional that ful- 
fills the constraints elaborated there will be very 
close to the energy functional one obtains from an 
Hamiltonian. 

• It should be mentioned that the use of the pro- 
jected density entering effective density-dependent 
vertices ^[p^] has already been proposed and used 
in Refs. [15-17]. However, in those references, a hy- 
brid approach is set up where transition densities 
are used in other parts of the functional, and for 
the restoration of spatial symmetries. 5 It has been 
pointed out in Ref. [13] that such hybrid approach 
may lead to unphysical results when set up for the 
restoration of spatial symmetries. Here, the theory 
is completely formulated in terms of the projected 
one- and two-body density matrices only. An open 
question that has to be addressed in the future is 
if and how the strategy to set up the functional 
we propose here can be generalized to the restora- 
tion of spatial symmetries, and perhaps even more 
general configuration mixing without becoming nu- 
merically intractable. 

IV. DISCUSSION AND CONCLUSION 

In this work, projection made by MR-EDF techniques 
including the recently proposed regularization [6, 8, 11] is 
further analyzed for the case of particle-number restora- 
tion of quasi-particle vacua of Bogoliubov type. Starting 
from a simple toy functional where the two-body effec- 
tive interaction is not explicitly density dependent, we 
show that the regularized energy can almost be regarded 
as a functional of the one- and two-body densities of the 
projected state. To follow the density functional method- 
ology, we propose to slightly modify the functional such 



5 For the special case of a pure particle-number projected MR-EDF 
calculation, the functional used in [17-19] could be mapped on 
a functional of the one- and two-body density matrices of the 
projected state. 



that it becomes a function of projected state densities 
only, and that the projected state becomes the inter- 
mediate trial state from which the functional and other 
observables are constructed. For particle-number projec- 
tion, such a modification could for instance be achieved 
within standard MR-EDF by slightly modifying the reg- 
ularization proposed in Ref. [6] while still removing 
the pathologies. Such an alternative interpretation may 
eventually serve as a justification of MR-EDF framework 
within a functional approach for particle number restora- 
tion when the effective kernels are not density dependent. 

As a matter of fact, most of the functional forms used 
nowadays do not enter into the class of functionals which 
can be safely used in MR-EDF. We show, however, that 
such functional can still be used in a symmetry restoring 
framework that does not make use of the MR-EDF tech- 
nique, but directly formulates the theory in terms of the 
one- and two-body density matrices of projected product 
states. 

This theory can be seen as a direct extension of the SR- 
EDF level that we proposed and is called here Symmetry- 
Conserving EDF approach. An illustration of the result- 
ing projected energy is given, showing that the method 
could be a valuable tool for the description of the ground 
state of a system with pairing including the restoration 
of particle number even when density dependence with 
non-integer powers is used in the functional. 

The analysis of similarities and differences between the 
MR-EDF theory and symmetry-conserving approaches 
was greatly simplified here because the original quasi- 
particle state and the projected state share the same 
canonical basis. For instance, expression (19) only holds 
in the canonical basis. In the present article, the ap- 
plications are restricted to projection after variation for 
which this equation is perfectly suited. The next the step 
will be the extension approach to perform variation after 
projection (VAP). VAP is usually solved using MR-EDF 
techniques by making variations with respect to the com- 
ponents of the original quasi-particle vacuum and not the 
projected state itself [20-23]. In the symmetry conserv- 
ing approach, one could follow the same strategy as in 
the standard MR-EDF approach, i.e. perform variations 
of the reference state. Work in that direction is currently 
underway. 

Last, we would like to mention that the present ar- 
ticle only discusses the case of particle-number projec- 
tion and the possibility to determine the ground-state 
energy. The MR-EDF technique is frequently used to 
restore other symmetries and to calculate excited states 
in a Generator-Coordinate framework. What these other 
configuration mixings have in common, is the fact that 
there does not exist a common canonical basis in which 
the one-body density matrices of the original and of the 
correlated states are simultaneously diagonal. An im- 
portant point to be clarified is if and how the formalism 
developed here can be generalized to those more general 
configuration mixings, and that without becoming nu- 
merically intractable. Finally, it has to be stressed that 
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the method proposed and explored here is not meant to functional is 
replace the MR EDF framework for the description of 
excited states and transition moments in complex nuclei. 
Instead, it might provide a numerically much more ef- 
ficient alternative to the MR EDF scheme when one is 
interested just in the ground state and its evolution, ei- 
ther in dynamics or thermodynamics. 
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where rq and k% refer to the occupation probabilities and 
anomalous densities of the state 4>o- Following Ref. [6], 
we decompose the energy kernels entering into the in- 
tegral of Eq. (cq:ekernel) into three terms £ p , £ pp and 
gnu corresponding to the kinetic, mean-field and pairing 
terms respectively. Then, £ pp and £ KK can be expressed 



v ijij n i n j 



£ pp = \Y: 

2 ^ 

ij 

+ ^E^>*M<5M<d , (A2) 

ij 

whereas 

CKK _ }_ -KK 0* 

ij 

+ 4 E "Stf ( K a 5 *jj [<P\ + K % dK a [<P]) 

ij 

+ 7E^J*«*M*««[V] ■ (A3) 



These expressions are the strict equivalent of the ones 
given in Eqns. (51-54) in Ref. [6]. For instance, regu- 
larizations have been obtained by removing terms with 
j = i in the last line of Eqs. (A2) and Eq. (A3). Accord- 
ingly, the spurious contribution to be removed from the 



& cg = \ E 1 'ijij J 5ni M Srii M ^n(0, <p)d<f , 

i 

£ CG = \ E »*« / 6 < M 5k * M W ' ^ ■ 



(A5) 

Therefore, when the regularization is included, this is 
equivalent to make the replacements 



Appendix A: Proof of Eq. (19) 

To prove Eq. (19), we have to explicitly remove terms 
that cause pathologies from the energy calculation as pro- an d 
posed in Ref. [6]. Starting from Eq. (50) of Ref. [6], the 
transition matrix elements can be expressed as 

n T = n °i + Sn i [<p] > 

K °a = K a + <*«« if] > 

= + 8^] , (Al) 



dtp n° v n° v J\f N (0,tp) 



2 71 



dtp (n ^ - Srii [tp] 5ni [ip] J Af N (0, ip) 



d( p { K Ti * K °ti - Sk *-i [<p] m) Aw(o, <p) 



in the last two terms of Eq. (13). 

From the equalities (Al), one can deduce new interest- 
ing relationships between the projected observables. For 
instance, performing the gauge integration of the first 
equation, we obtain 



= / d<pn o i ' p AfN(0,<p) = n i +Sni 



with 



(A6) 



Sn z = < - nj = / dtp 5m [tp] Af N (0, tp) . (A7) 



From this, let us now re-express the different quantities 
entering in Eq. (13) 



dtp n° j ' p jV N (0,tp)=n i n^ 



13 
2tv 



+ / dtp Srii [tp] 6rij[(p]J\fN{0,(p) 



where, in the specific case i = j, we recognize the term 
that enters in the regularization to be the last one. There- 
fore, the term entering into the regularization of £ pp can 
be expressed as 



dtp (n ^ n ^ - Srii [tp] Srii [tp]) Ajv(0, tp) 

= n°n° + 2n° l 8n l 

= nf nf ~ SmSm . (A8) 
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To derive an expression of the term entering in the 
regularization of £ KK , one can proceed in a similar way. 
We first define Sk*^ and Sku through 



dp <°*A6v(0, tp) = dtp (4* + 5k% [tp]) M N (0, if) 
Jo 

= K a + 3 K it 

T />27T 

dip K ^Af N (0, tp) = d v (4 + Sks [<p]) AT n (0, p) 



Therefore the term entering in the regularized functional 
is given by 

d< P {^m K °m ~ 5k % if] 5k s if]) N N (0, tp) 



One can then take advantage of the fact that 



(A9) 



2 71 



dip Srii [tp] Srij [ip]Af N (0,ip) 



I 0* , r * , 0* c 
t K s f% + OKtfKfo + K a OKii 

d<p8n% [ip] 5 [f] A/iv(0, ip) 



and that 

Srii [tp] Srii [tp] = —Sku [ip] <5k*j [ip] . (A10) 

The first equality is nothing but Eq. (18), whereas the 
second equality can be proved by expressing Srii [</?], 
5kh [ip] and 5k* s [tp] directly in terms of the Ui and Vi 
of the SR-EDF theory and the gauge angle tp starting 
from Eq. (72-74) of Rcf. [6]. Altogether, we obtain: 

( k«° K °it ~ 5k b [<P\ Sk » M ) Mv(0, tp) = nf 



dtp [r£ v - Srii [tp] Srii [tp] ) 7Vjv(0, tp) 



= (n?(l-n?) + 6m5m). (All) 

Combining this expression with Eq. (A8), we finally de- 
duce the expression (19) for the regularized functional. 



Appendix B: Particle number restoration with 
recurrence relation 



A method, alternative to the the gauge-integration 
method is presented here to calculate the one- and two- 
body density matrix components of a projected product 
state. This method turns out to be very fast and efficient 
numerically. 

Let us start from a quasi-particle state written in its 
canonical basis as 



$ °> = nU + ^44) |o>, 



(Bl) 



i>0 



where \xi\ 2 = n°/(l — n°). The associated projected state 
with N particles can be expressed as 

N 

(B2) 



\i>0 J 



Starting from these expressions, it has recently been 
shown [14] that the elements of the one- and two-body 
density matrix are given by 



n? = N\x^ lN - l{j) 



In 



'JV-l 



(B3) 

for (i^j), (B4) 
lN-i{i,j) 



IN 



(B5) 



while as already mentioned Rf m = nf . The different 
coefficients entering in n N and R N are given by: 



Ik = 

Ik® = 
lK(i,j) = 



i K )jti \ x h\ 2 



\ Xi i I 2 ' ' ' \ x iK I 



Direct use of these expressions for K = N is rather diffi- 
cult numerically. However, these coefficients verify sim- 
ple recurrence relations that are straightforward to im- 
plement on a computer. These recurrence relations have 
been recently used to solve numerically the Variation Af- 
ter Projection (VAP) [24, 25] and to set up a new func- 
tional for pairing accounting for particle-number conser- 
vation [14]. 

In the present work, we use the recurrence method 
to perform PAV within the symmetry-conserving EDF 
framework. In that case, a preliminary SR-EDF calcula- 
tion is performed leading to a quasi-particle state given 
by (Bl) with a set of {xi} values. Here, we have used the 
ev8 code [26]. From this set, the quantities In-i(i) and 
In arc evaluated via the recurrence relations 



IkH) 
Ik 



-(K-l)] Xi \ 2 I K -i(i) 
Xi \ 2 I K -i-(K -2) 



\x t \ 4 I K ~2{i) 



(B6) 



with the condition 7q = Io(i) = 1, -^l = ^2k\ x k\ 2 an d 



h(i) = h 



The occupation numbers of the pro- 



jected state can then be calculated as well as the corre- 
lation components using the relation [14, 27]: 



I^iljj x i ■''.! 



nf ~ nf 



ijij 



for (B7) 
(B8) 

where for i = ?', we have = and R^n — 0. This 
method is referred to as "recurrence method" in the text. 
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